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D . Abstract 

Energy-efficient and flexible communication protocols are instrumental for effective and pervasive 
development of control and monitoring applications over wireless networks. The IEEE 802.15.4 protocol 
stack is the de-facto reference communication standard for industrial and home automation. Although 
the performance of the medium access control (MAC) of this protocol has been thoroughly investigated 
under the assumption of ideal wireless channel, there is still a lack of understanding of the cross- 
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layer dynamics between MAC and physical layer in the presence of realistic wireless channel models 



that include path loss, multi-path fading and shadowing. An investigation of these mutual effects is 
instrumental for consistent performance prediction, correct design and optimization of the protocols in 
■ the prospected application scenarios. In this paper, a novel approach to analytical modeling of these 

o 

_ interactions and performance prediction is proposed. The analysis considers simultaneously a composite 



channel fading, interference generated by multiple terminals, the effects induced by hidden terminals, and 
the MAC reduced carrier sensing capabilities. New results on the MAC-physical layer interactions are 
derived and validated for single-hop and multi-hop topologies. It is shown that performance indicators 
of the MAC over fading channels are often far from those derived under ideal channel assumptions. 
Moreover, it is established to what extent fading may be beneficial for the overall network performance 
and how the results presented in this paper can be used to drive joint optimization of physical layer 
and MAC layer parameters. 
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I. Introduction 

The development of wireless sensor network (WSN) systems relies heavily on understanding 
the behavior of underlying communication mechanisms. When sensors and actuators are inte- 
grated within the physical world with large-scale and dense deployments, potential mobility of 
nodes, obstructions to propagation, fading of the wireless channel and multi-hop networking 
must be carefully addressed to offer reliable services. For instance, this is essential for wireless 
control applications. Indeed, wireless interfaces can represent bottlenecks as they may not provide 
links as solid as required by applications in terms of guarantees on reliability, delay, and energy 
efficiency. 

There seems to be consensus on that the protocols for physical layer and medium access 
control (MAC) for low data rate and low power applications in the future will be based on the 
flexible IEEE 802.15.4 standard with its numerous variants [1J. That standard has been indeed 
adopted with some modification also by a number of other protocol stacks, including ZigBee, 
WirelessHART, ISA-100. It is already being used for applications in industrial control, home 
automation, health care, and smart grids. Nevertheless, there is not yet a clear understanding 
of the achievable performance of the IEEE 802.15.4 protocol stack, with the consequent in- 
ability to adapt the communication performance (e.g., through cross-layer optimization) to meet 
challenging quality of service constraints. 

The IEEE 802.15.4 MAC layer has received much attention, with focus on performance 
characterization in terms of reliability (i.e., successful packet reception probability), packet delay, 
throughput, and energy consumption. Some initial works, such as 10, are based on Monte Carlo 
simulations. More recent investigations have attempted to model the protocol performance by 
theoretical analysis for single hop networks [[3]|-[|9]|. These analytical studies are almost all 
based on extensions of the Markov chain model originally proposed by Bianchi for the IEEE 
802.11 MAC protocol IfTOl . In [fTTI . modelling of packet losses due to channel fading have been 
introduced into the homogeneous Markov chain developed for the IEEE 802.15.4 MAC setup 
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presented in [5]. However, fading is considered only for single packet transmission attempts, 
the effect of contention and multiple access interference is neglected, and the analysis is neither 
validated by simulations nor by experiments. In lfT2l . [fT3l . the optimal carrier sensing range is 
derived to maximize the throughput for IEEE 802.11 networks; however, statistical modeling of 
wireless fading has not been considered, but a two-ray ground radio propagation model is used 
in [fT2l and only the path loss is considered in [Tl3l . 

The technical contributions surveyed above focus on modeling the typical MAC mechanisms, 
such as the finite retry limits of the carrier sense multiple access with collision avoidance 
(CSMA/CA), the carrier sensing range, saturated and unsaturated traffic conditions. However, as 
already stated, the presence of the wireless channel has been either neglected or only marginally 
considered. Actually, simulations show that the studies mentioned above may yield highly 
inaccurate results in terms of reliability, delay, and energy consumption as soon as realistic 
wireless channel behaviours and interference profiles are introduced. In this regard, recent studies 
have investigated the performance of multiple access networks in terms of multiple access 
interference and capture effect lfT4l . However, in [14j the model of the MAC mechanism is limited 
to homogeneous single-hop networks (same wireless channel for every node) with uniform 
random deployment. Moreover, the fading caused by multi-path propagation and shadowing 
effects have not been considered, while it is known that it may have a crucial impact on the 
performance of packet access mechanisms [fT5l . 

In this paper we propose a novel analytical model that is able to capture the cross-layer inter- 
actions of IEEE 802.15.4 MAC and physical layer over interference-limited wireless channels 
with composite fading models. The main original contributions are as follows: 

• We propose a general modeling approach for characterization of the MAC performance 
with heterogeneous network conditions, a composite Nakagami-lognormal channel, explicit 
interference behaviours and cross-layer interactions. 

• Based on the new model, we determine the impact of fading conditions on the MAC 
performance under various settings for traffic, inter-node distances, carrier sensing range, 
and signal-to-(interference plus noise)-ratio (SINR). We show how existing models of the 
MAC available from the literature may give unsatisfactory or inadequate predictions for 
performance indicators in fading channels. 
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• We discuss system configurations in which a certain severity of the fading may be beneficial 
for overall network performance. Based on the new model, it is then possible to derive 
optimization guidelines for the overall network performance, by leveraging on the MAC- 
physical layer interactions. 

To determine the network operating point and the performance indicators in terms of reliability, 
delay, and energy consumption for single-hop and multi-hop topologies, a moment matching 
approximation for the linear combination of lognormal random variables based on Ifl6l and lfT71l 
is adopted in order to build a Markov chain model of the MAC mechanism that embeds the 
physical layer behavior. In particular, the challenging part of the new analytical setup proposed 
in this paper is to properly model the rather complex interaction between the MAC protocol 
and the wireless channel with explicit description of the dependence on several topological 
parameters and network dynamics. For example, we include failures of the channel sensing 
mechanism and the presence of hidden terminals, namely nodes that are in the communication 
range of the destination but cannot be listened by the transmitter. Whether two wireless nodes 
can communicate with each other depends on their relative distance, the transmission power, 
the wireless propagation characteristics and interference caused by concurrent transmissions on 
the same radio channel: the higher the SINR is, the higher the probability that packets can be 
successfully received. The number of concurrent transmission depends on the traffic and the 
MAC parameters. In this paper, we describe how to account for statistical fluctuations of the 
SINR in the Markov chain model of the MAC. 

The remainder of the paper is organized as follows. In Section HH we introduce the network 
model. In SectionUm we derive an analytical model of IEEE 802.15.4 MAC over fading channels. 
In Section [TV] reliability, delay, and energy consumption are derived and the results are extended 
to multi-hop topologies in Section [V] The accuracy of the model is evaluated by Monte Carlo 
simulations in Section [VH along with a detailed analysis of performance indexes with various 
parameter settings. Section IVIII concludes the paper and prospects our future work. 

II. Network Model 

We illustrate the network model by considering the three topologies sketched in Fig. [Q 
Nevertheless, the analytical results that we derive in this paper are applicable to any fixed 
topology. 
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The topology in Fig. [Tk) refers to a single-hop (star) network, where node i is deployed at 
distance r i from the root node at the center, and where nodes forward their packets with single- 
hop communication to the root node. The topology in Fig. []})) is a multi-hop linear topology, 
where every node generates and forwards traffic to the root node by multi-hop communication. 
The distance between two adjacent nodes is denoted as r it j. In Fig. [it), we illustrate a multi- 
hop topology with multiple end-devices that generate and forward traffic according to an uplink 
routing policy to the root node. 

Consider node i that is transmitting a packet with transmission power P tx ,i- We consider an 
inverse power model of the link gain, and include shadowing and multi-path fading as well. The 
received power at node j, which is located at a distance r it j, is then expressed as follows 

P™,i,3 = C ° P ^ 1 fi exp(yi) . (1) 

T i,j 

The constant c represents the power gain at the reference distance 1 m, and it can account for 
specific propagation environments and parameters, e.g., carrier frequency and antennas. In the 
operating conditions for IEEE 802.15.4 networks, the inverse of c (i.e., the path loss at the 
reference distance) is in the range 40 — 60 dB 0]|. The exponent k is called path loss exponent, 
and varies according to the propagation environment in the range 2 — 4. The factor fi models 
the channel fading due to multi-path propagation, which we assume to follow a Nakagami 
distribution with parameter k > 0.5 and p.d.f. 

where T(k) is the standard Gamma function T(k) = exp(— x)x K ~ 1 dx. We remark here 
that the Nakagami distribution is quite general and can model fading environment with various 
degrees of severity. A random lognormal component models the shadowing effects due to 
obstacles, with A/"(0,of). The standard deviation <jj is called spread factor of the shadowing. 

The model we adopt is accurate for IEEE 802.15.4 in a home or urban environment, in which 
devices may not be in visibility. In the following section, a generalized model of a heterogeneous 
network using unslotted IEEE 802.15.4 MAC over multi-path fading channels is proposed. 
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III. IEEE 802.15.4 MAC and PHY Layer Model 

In this section we propose a novel analytical setup to derive the network performance indica- 
tors, namely the reliability as probability of successful packet reception, the delay for successfully 
received packets, and the average node energy consumption. We first consider a single-hop case, 
and then we generalize the model to the multi-hop case in Section |V] 

A. Unslotted IEEE 802.15.4 MAC Mechanism 

According to the IEEE 802.15.4 MAC, each node can be in one of the following states: (i) 
idle state, when the node is waiting for the next packet to be generated; (ii) backoff state; (iii) 
clear channel assessment (CCA) state; (iv) transmission state. 

Let a node be in idle state with probability pidi e ,i- The three variables given by the number 
of backoffs NB, backoff exponent BE, and retransmission attempts RT are initialized: the 
default initialization is NB=0, BE=m , and RT=0. From idle state, a node wakes up with 
probability which represents the packet generation probability in each time unit of duration 
Sb = aUnitBackoffPeriod, and moves to the first backoff state, where the node waits for a random 
number of complete backoff periods in the range [0, 2 BE — 1] time units. 

When the backoff period counter reaches zero, the node performs the CCA procedure. If the 
CCA fails due to busy channel, the value of both NB and BE is increased by one. Once 
BE reaches its maximum value m b , it remains at the same value until it is reset. If NB 
exceeds its threshold m, the packet is discarded due to channel access failure. Otherwise the 
CSMA/CA algorithm generates again a random number of complete backoff periods and repeats 
the procedure. A node is in CCA state with probability r-i, and either moves to the next backoff 
state if the channel is sensed busy with probability or moves to transmission state with 
probability (1 — c^). The node experiences a delay of aTurnaroundTime to turn around from 
listening to transmitting mode. 

The reception of the corresponding ACK is interpreted as successful packet transmission. The 
node moves from the transmission state to idle state with probability (1 — 7ij). As an alternative, 
with probability 7^, the packet sent to node j is lost, and the variable RT is increased by one. 
As long as RT is less than its threshold n, the MAC layer initializes BE=m and starts again 
the CSMA/CA mechanism to re-access the channel. Otherwise the packet is discarded due to 
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the retry limit. 



B. MAC-Physical Layer Model 

In this subsection, we substantially extend the MAC model presented in lfT8l . which was 
developed under ideal channel conditions, to include the main features of channel impairments 
and interference. 

Assume packets are generated with Poisson distribution at rate A/. The probability of generation 
of a new packet after an idle unit time is then qi = 1 — exp(— A//S&). Effects of a limited buffer 
size can be included in the expression of qi, by considering the probability that the node queue 
is not empty after a packet has been successfully sent, after a packet has been discarded due to 
channel access failure or due to the retry limit (see [fT8l ). 

We define the packet successful transmission time L s and the packet collision time L c as 



L s — L + t ack + L ack + IFS . 

Lc L -\- t m ac k j 



(2) 



where L is the total length of a packet including overhead and payload, t ack is ACK waiting 
time, L ack is the length of ACK frame, IFS is the inter-frame spacing, and t m>ack is the timeout 
(waiting for the ACK) in the retransmission algorithm, see details in [0Q. 

The expression of the idle state probability p\d\ c ,i can be derived from the normalization 
condition of the Markov chain model in |fT8l , as 
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where yi = 7ij(l — ctT +1 )- The CCA probability t; can be derived as 



The busy channel probability is 

Oti = a pk t,j + «ack,i , (5) 

where a pktyi is the probability that node i senses the channel and finds it occupied by an ongoing 
packet transmission, whereas a ac k,j is the probability of finding the channel busy due to ACK 
transmission. We are now ready to include the effects of fading channel as it follows. 

The busy channel probability due to packet transmissions evaluated at node % is the combination 
of three events: 

• at least one other node has accessed the channel within one of the previous L units of time; 

• at least one of the nodes that accessed the channel found it idle and started a transmission; 

• the total received power at node i is larger than a threshold a, so that an ongoing transmission 
is detected by node i. 

The combinations of all busy channel events yields the busy channel probability that node i 
senses the channel and finds it occupied by an ongoing packet transmission 

a^ = LU,{P? ct ) , (6) 



where 
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and 



pdet = Pr 



m+1 



(8) 



is the detection probability. The index k accounts for the events of simultaneous accesses to the 
channel and the index j enumerates the combinations of events in which a number I of channel 



accesses are performed in the network simultaneously. Given N nodes in the network, the index 
kj refers to the node in the k-th position in the j-th combination of / out of N — 1 elements 
(node i is not included). The index q n accounts for the combinations of events in which one or 
more nodes find the channel idle simultaneously. 

The busy channel probability due to an ACK transmission follows from a similar derivation. 
An ACK is sent only after a successful packet transmission: 



ack, l = ^ack^((l- 7(/ „, 1 ,)^ dCt ) , (9) 



a 



where L ack is the length of the ACK. The index w denotes the destination node of q n in the 
expression of Hi. The detection probability Pf et is evaluated with respect to the set of destination 
nodes. By summing up Eqs. © and ©, we compute a* in Eq. ©. 

We next derive an expression for the packet loss probability 7^, namely the probability that a 
transmitted packet from node i is not correctly detected in reception by node j, by using similar 
arguments as above. A packet transmission is not detected in reception if there is at least one 
interfering node that starts the transmission at the same time and the SINR between the received 
power from the transmitter and the total interfering power plus the noise level N is lower than a 
threshold b (outage). In the event of no active interferers, which occurs with probability l-H^l), 
the packet loss probability is the probability that the signal-to-noise ratio (SNR) between the 
received power and the noise level is lower than b. Hence, 

7 y = (1 - 7U{1)) Ptf + % (Pif) + (2 L - 1) % ((1 - P, dct ) P°?) , (10) 

where P^ d is the outage probability due to composite channel fading on the useful link (with 
no interferers), 
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and is the outage probability in the presence of interferers (with composite channel fading 
on every link), 



pout _ p 
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(12) 
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The expressions of the carrier sensing probability r\ in Eq. ©, the busy channel probability 
oti in Eq. ©, the collision probability in Eq. (flOl) , for I — 1, . . . , N, form a system of non-linear 
equations that can be solved through numerical methods |fl9ll . 

We next need to derive the detection probability and the outage probabilities in the devised 
wireless context. With such a goal in mind, we present some useful intermediate results in the 
following section. 

C. Model of Aggregate Multi-path Shadowed Signals 

In this section, we approach the problem of computing the sum of multi-path shadowed signals 
that appear in the detection probability and in the outage probability. 

Consider node i performing a CCA and let us focus our attention on the detection probability in 
transmission Pr En=i ^rx,n,t > a ]> where x is the current number of active nodes in transmission. 
By recalling the power channel model in Eq. (Q~|), let us define the random variable 

Yi = In I Y]A i>n exp(y ri 

\n=l 

with 

■™-i,n Co ^ Jni 

and y n ~ A/"(0, cr^). Since a closed form expression of the probability distribution function of 
Yi does not exist, we resort to a useful approximation instead. In order to characterize Yi, we 
apply the Moment Matching Approximation (MMA) method, which approximates the statistics 
of linear combination of lognormal components with a lognormal random variable, such that 

Yi~Af( m ,a 2 Y J. 

According to the MMA method, r] Yi and a Yi can be obtained by matching the first two moments 
of exp(Yj) with the first two moments of J2n=i A,nexp(i/ n ), i.e., 

M x = exp y-r) Yi + ^(T Y ^j = E{Ai, n } exp (^ yn + ^a yn j , (13) 



M 2 4 exp (-2% + 2(7*,) =Y1 E {AmAi, n } expL ym -H] yn +(^ + ^f+ Py^y-A 

m=l n=l ^ ^ 

(14) 
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Solving Eqs. (Q]§ and (HU) for r/y. and oy. yields r/y. = 0.51n(M 2 ) - 21n(Mi), and erf, 
ln(M 2 ) -21n(M x ). It follows that 

'In (a) - 7/yj 



jdet 



Pr 



/ - 1 rx,n,'. 



> a 



n=l 



where 



Pr [exp(Fj) > a] « Q 



(15) 



Similar derivations follow for the outage probability in reception: 
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which yields rjy . = 0.5 ln(M 2 ) - 2 lnfMA cr| = ln(M 2 ) - 2 ln(Mi). Therefore, 
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The analysis above holds for a generic weighted composition of lognormal fadings. In the 
case of lognormal channel model, where only shadow fading components are considered, (i.e., 
fi = 1), the outage probability becomes 



pout _ p 
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ln(fe) — fly 
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For a Nakagami-lognormal channel, the outage probability becomes 
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For integer values of m, the integration in z yields 

Jo V2x(T 9iij w \ 2°y. J ^ T{i + 1) 

The mean and standard deviation of Yi and Y^j can be obtained by inserting the moments of /, 
in the moments of A i>n and -Bjj>. For Gamma distributed components fi, we obtain E{fi} = 1 
and E{f?} = a±l. 

We remark here that the evaluation of Pf et , and P°J lt can be carried out off-line, with respect 
to the solution of the system of nonlinear equations that need to be solved when deriving 77, 
a,i and 7j j. Therefore, the proposed model can be implemented with only a slight increase of 
complexity with respect to the analytical model of the IEEE 802.15.4 MAC mechanism presented 
in [fT8l , and the online computation time is not affected significantly. 



IV. Performance Indicators 

In this section, we investigate three major indicators to analyze the performance of the 
IEEE 802.15.4 MAC over fading channels. These indicators will also be used to validate the 
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analytical model we derived in the previous section, by comparing results obtained from the 
(approximate) model with those obtained by extensive simulation campaigns. The first one is 
the reliability, evaluated as successful packet reception rate. Then we consider the delay for 
the successfully received packets as the time interval from the instant the packet is ready to 
be transmitted, until an ACK for such a packet is received. Eventually, we consider the energy 
consumption of network nodes. 

A. Reliability 

For each node of the network, the reliability is based on the probability that packets are 
discarded at MAC layer. In unslotted CSMA/CA, packets are discarded due to either (i) channel 
access failure or (ii) retry limits. A channel access failure happens when a packet fails to obtain 
clear channel within m + 1 backoff stages. Furthermore, a packet is discarded if the transmission 
fails due to repeated packet losses after n + 1 attempts. Therefore, the reliability is 

Ri = 1 — P c f,i — Pa,i • (18) 

According to the IEEE 802.15.4 MAC mechanism described in Section Hll-A[ the probability 
that the packet from node i to node j is discarded due to channel access failure can be expressed 
as 

p «r +1 (i-(7M-(i-< +1 )) n+1 ) 
i- 7M (i-«r +1 ) 

and the probability of a packet discarded due to retry limits is 

p cr , = (7M(i-«r +1 )) n+1 - 

It is worthwhile mentioning that the last expressions link the reliability at the MAC level with 
the statistical description of wireless channel environment through Eq. (flOl) and the analysis of 
Subsection IIII-CI 



B. Delay 

We define the delay Di for successfully delivered packets from node i. If a packet is discarded 
due to either the limited number of backoff stages m or the finite retry limit n, its delay is not 
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included into the average delay. 

Let D ijh be the delay for node i that sends a packet successfully at the h-th attempt. The 
expected value of the delay is 

n 

E[A]= EPrFWE[A,ft], ^ 19 ) 

/i=0 

where the event % t denotes the occurrence of a successful packet transmission at time h + 1 
given /i previous unsuccessful transmissions, whereas the event indicates a successful packet 
transmission within n attempts. Therefore, we can derive 

Pr \% t m = ^ 1 j- = ± ^ - '-L^ — — - L . (20) 

£Lo(7m(i-< +1 ))* i-(7M(i-«r +1 )) n+1 

We recall that is the collision probability, which we derived in Eq. (flOl) together with 
Eqs. (TT3T) and (fTTI) . and 1 — a 4 m+1 is the probability of successful channel access within the 
maximum number of m backoff stages, where ct™ +1 follows from ©. 
The average delay at the h-th. attempt is 

h 

E[D lth ] = L s + hL c + J2 E lTi], (21) 

z=o 

where TJ is the backoff stage delay, whereas L s and L c are the time periods in number of time 
units for successful packet transmission and collided packet transmission computed in Eq. (0. 

Since the backoff time in each stage k is uniformly distributed in [0, Wk — 1], where Wk = 2 BE , 
the expected total backoff delay is 

m ( r Wh - 1 \ 

E[T h ] =T SC + J2 (rT sc + J2^—S b , (22) 

r=0 V A:=0 ^ / 

where T sc is the sensing time in the unslotted mechanism. The event 3t r denotes the occurrence 
of a busy channel for r consecutive times, and then an idle channel at the (r + l)th time. By 
considering all the possibilities of busy channel during two CCAs, the probability of @> r is 
conditioned on the successful sensing event within m attempts given that the node senses an 
idle channel in CCA. It follows that 
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By applying Eqs. (I2TT) - (|23l ) in Eq. (fT9l) , the average delay for successfully received packets 
is computed. Note that the delay is experienced at the MAC level and is hereby linked to the 
fading channel through the dependency on oti and 7,^- evaluated in the previous section. 

C. Energy Consumption 

Here we develop the expression of the total energy consumption of node i as the sum of the 
contributions in backoff, carrier sense, transmission, reception and idle-queue state: 

Etot,i = E bii + E 8i i + E t)i + E r> i + Eq }i . (24) 

In the following, each component of this formula is derived according to the state probabilities 
in Section IIII-AI The energy consumption during backoff is 

(l-(2aQ m+1 )(l-a i ) \ 
(l-2^)(l-«r +1 ) J ' 

where P id is the average power consumption in idle-listening state, as we assume that the radio is 
set in idle-listening state during the backoff stages. The energy consumption for carrier sensing 
is E Sji = P sc Ti, where P sc is the average node power consumption in carrier sensing state. The 
energy consumption during the transmission stage, including ACK reception, is 

E t ,i =(1 - a^niPtL + P id + (P r (l - 7i)i ) + Pid7ij)^ack) , 

where P t and P r are the average node power consumption in transmission and reception respec- 
tively, and we assume t miac k = -^ack + 1 in backoff time units S&. In the single-hop case, we 
assume that the node is in sleeping state with negligible energy consumption during inactivity 
periods without packet generation. Hence, the energy consumption during the idle-queue state 
is given by 

Et,i =P s Pid\e,i , 

where P s is the average node power consumption in sleeping mode, and Pidi e ,i is the stationary 
probability of the idle-queue state as derived in Eq. ©. 

In the multi-hop case, as we will see in the next section, relay nodes are in idle-listening 
state also during the inactivity period (because of the duty cycle policy), and an extra cost for 
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receiving packets and sent ACKs has to be accounted for. 

V. Multi-hop Communications 

Here we extend the Markov chain model to a general network in which information is 
forwarded through a multi-hop communication to a sink node. 

The model equations derived in Section IIII-BI are solved for each link of the network, by 
considering that a generic node i forwards aggregate traffic Qi. The total average aggregated 
traffic of node i is Qi = qi/S b pkt/s, where is the probability of having a packet to transmit 
in each time unit and S& is the duration of the basic time unit in IEEE 802.15.4. The aim of the 
following analysis is to investigate the total traffic load Q h which we associate to the probability 
qi in the per- link model equations. To do so, we need to characterize the traffic distribution in 
the network according to a routing policy. 

We define A = [0, Ai, \n] the vector of node traffic generation rates, where each component 
is associated to a node. In addition to Ai, node i has to forward traffic generated by its children 
nodes. The effect of routing can be described by the routing matrix M G R( iV + 1 ) x ( iV + 1 ), in 
which Mjj = 1 if node j is the destination of node i, and M^j = otherwise. We assume that 
the routing matrix is built such that no cycles exists. We define the traffic distribution matrix T 
by scaling M by the probability of successful reception in each link as only successfully received 
packets are forwarded, i.e., 

I i. j H,j Ri j 

where Ri is given by Eq. (fT8l . 

Therefore, the vector of node traffic generation probabilities Q = [0,Qi, ■ ■ ■ ,Qn] is the 
solution of a system of flow balance equations Q = Q T + A. In steady state, we have 

Q = A[I-T]- 1 . (25) 

where Ig R( 7V+1 ) x ( Ar+1 ) is the identity matrix. As routing is acyclic, it is easy to show that T 
has spectral radius less than one, therefore the inverse always exists. 

Eq. (|25T) gives the relation between the idle packet generation probability q if the routing matrix 
M, and the performance at MAC layer (through the link reliability Ri). To obtain the multi- 
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hop network model, we couple Eq. (|25T ) with the expressions for r», a, and Ri, as obtained by 
Eqs. ©, ©, and 4TB). 

With respect to single-hop networks, there is an extra energy consumption E x>i due to the 
packets and ACKs received by relay nodes based on the routing matrix M, 

N 

E x ,i = 22 M nA l ~ 7n,i)(l - a n )r n (P t L + P id + (P r (l - y n j) + P\dln,i)L^ k ) . 

n=l 

We validate and exploit these analytical results in the next section. 

VI. Performance Evaluations 

In this section, we present numerical results for the new model for various settings, network 
topology, and operations. We report extensive Monte Carlo simulations to validate the accuracy 
of the approximations that we have introduced in the model. The settings are based on the 
default specifications of the IEEE 802.15.4 [1J. We perform simulations both for single-hop and 
multi-hop topologies. 

A. Single-hop Topologies 

In this set of performance results, we consider a single-hop star topology in Fig. [T^). We let 
the number of nodes be N = 7, the MAC parameters m = 3, m = 4, m b = 7, n = 0, L = 7, 
L^k = 2 and the physical layer parameters P tx ,i = dBm, and k = 2. We validate our model 
and study the performance of the network by varying the traffic rate Aj = A, % = 1, N, in the 
range 0.1 — 10 pkt/s, the radius r i = r, i — 1, ...,N, in the range 0.1 — 10 m, the spread of 
the shadow fading cr, = a, i = 1, N, in the range — 6, and the Nakagami parameter k in 
the range 1 — 3. Moreover, we show results for different values of the carrier sensing threshold 
a = —76, 66, 56 dBm, and outage threshold b = 6, 10, 14 dB. 

In Fig. |2l we report the average reliability over all links by varying the node traffic rate A. 
The results are shown for different values of the spread a in the absence of multi-path (fa = 1). 
There is a good matching between the simulations and the analytical model (TTSl . The reliability 
decreases as the traffic increases. Indeed, an increase of the traffic generates an increase of the 
contention level at MAC layer. However, we can observe that the impact of shadow fading can 
be more relevant with respect to variations in the traffic. Therefore, a prediction based only on 
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Markov chain analysis of the MAC without including the channel behavior, as in the previous 
literature, is typically inaccurate to capture the performance of IEEE 802.15.4 wireless networks, 
especially at larger shadowing spreads. 

In Fig. [3] the average delay over all links is reported. Also in this case simulation results follow 
quite well results obtained from the model as given in Eq. (fT9l ). An increase of traffic leads to an 
increase of the average delay due to the larger number of channel contentions and consequently 
an increase in the number of backoffs. The spread of shadowing components does not impact 
on the delay significantly, particularly for low traffic, because lost packets due to fading are 
not accounted for in the delay computation. When the traffic increases, we note that fading is 
actually beneficial for the delay. In fact, the delay of successfully received packets reduces by 
increasing a. This is because the occurrence of a deep fading reduces the probability of successful 
transmission. However, since this holds for all nodes, the average number of contending nodes 
for the CCA may reduce, thus reducing the average delay of successfully received packets. It 
is not possible to capture this network behavior by using separate models of the IEEE 802.15.4 
MAC and physical layers as in the previous literature, since this effect depends on a cross-layer 
interaction. 

In Fig. HI the average power consumption over all links is presented and compared to the 
analytical expression in Eq. (1241) . The number of packet transmissions and ACK receptions is 
the major source of energy expenditure in the network. Therefore, an increase of the traffic 
leads to an increase of the power consumption, while performance are marginally affected by 
the spread of the fading. However, the power consumption is slightly reduced when the spread 
is a = 6, due to the smaller number of received ACKs. Note that no power control policy is 
implemented. 

In Fig. [5l the average reliability is reported as a function of the radius r for different values 
of the spread o. Again, analytical results, obtained through Eq. (fT8~l) . are in good agreement with 
those provided by simulations. For the ideal channel case (i.e., a = 0) the size of the network 
does not affect the reliability in the range r = 0.1 — 10 m. For a = 6, the performance degrades 
significantly as the radius increases. An intermediate behavior is obtained for a = 3, where the 
reliability is comparable to the ideal channel case for short links, but it reduces drastically for 
r > 1 m. The effect is the combination of an increase of the outage probability with the radius 
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(due to the path loss component), and hidden terminals that are not detected by the CCA. 

In Fig. [6l we report the average delay by varying the radius r for different values of the 
spread o. The shadowing affects the delay positively and the effect is more significant for 
larger inter-node distances: in this case the average number of contending nodes for the free 
channel assessment reduces, thus the busy channel probability reduces, which in turn decreases 
the average delay of successfully received packets. 

In Fig. [7J the average power consumption by varying r is presented. We notice a similar 
behavior as for the delay. The power consumption reduces with the fading and the increasing 
size of the network. Nodes spend less time in the backoff and channel sensing procedure due to 
reduced number of contending nodes and the number of ACKs. 

Fig. [8] shows the average reliability as a function of the shadowing spread a. The results are 
plotted for different values of the carrier sensing threshold a. The reliability decreases when 
the threshold a become larger. The impact of the variation of the threshold a is maximum 
for a = 0, and the gap reduces when the spread a increases. In Fig. [9l the average delay 
is plotted as a function of the spread a. Depending on the threshold a, the delay shows a 
different behavior when increasing a: it increases for a = — 76 dBm and it decreases for a = 
—66 dBm, and a = —56 dBm. As we discussed above, the spread a may reduce the delay under 
some circumstances. However, when the threshold is large, the average number of contenders 
is less influenced by the fading and does not decrease significantly, while the busy channel 
probability becomes dominant and the number of backoffs increases, so that the delay increases 
as well. Fig. \T0\ reports the average power consumption by varying the spread a. The power 
consumption reduces by increasing the threshold a as a consequence of the smaller number of 
ACK transmissions, although a maximum consumption is observed for low values of the spread 
a. 

In Fig. [HI we plot the average reliability as a function of the spread a for different values of 
the outage threshold b. The threshold b does not affect the performance noticeably for o = 0, 
while the gap in the reliability increases with a. Note that for a high threshold the reliability 
tends to increase with a as long as a is small or moderate, and it decreases for large spreads. 
In our setup, a maximum in the reliability is obtained for a w 2. 

In Fig. [121 we report the combined effects of shadow fading and multi-path fading on the 
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reliability. We show the reliability as a function of the spread a of the shadow fading for different 
values of the Nakagami parameter k. We recall that k — 1 corresponds to Rayleigh fading. There 
is a good match between the simulations and the analytical model (fl~8l) . The effect of the multi- 
path is a further degradation of the reliability. However, the impact reduces as the Nakagami 
parameter k increases and the fading becomes less severe. In fact, for k 1, the effect of 
multi-path becomes negligible. Furthermore, the multi-path fading and the composite channel 
evidences the presence of the maximum at a w 2 in the plot of reliability. 

B. Multi-hop Linear Topologies 

In this set of performance results, we consider the multi-hop linear topology in Fig. The 
number of nodes (and hops) is iV = 5, with the same MAC and physical layer parameters as 
in the single-hop case. We validate our model and study the performance of the network as a 
function of the hop distance r^- in the range r = 0.1 — 10 m, and the spread of the shadow 
fading in the range a = — 6. We show results for each hop, and for different values of the 
carrier sensing threshold a = —76, 66, 56 dBm, and outage threshold 6 = 6, 10, 14 dB. 

In Fig. \T3\ the end-to-end reliability is reported from each node to the destination node for 
different values of the spread a. The analytical model follows well the simulation results. The 
end-to-end reliability decreases with the hop number. This effect is more evident in the presence 
of shadowing. Fig. [14] shows the end-to-end reliability from the farthest node to the destination 
by varying the distance r between every two adjacent nodes for different values of the spread a. 
The reliability is very sensitive to an increase of hop distance. In Fig. [T31 we show the end-to-end 
reliability by varying the spread a of the shadow fading. Results are shown for different values 
of the carrier sensing threshold a. In Fig. [161 we plot the end-to-end reliability for different 
values of b. Similar considerations as for the single-hop case applies here. However, for the 
linear topology, the reduction of the carrier sensing range from a = — 76 dBm to a = —66 dBm 
influences less the reliability since hidden nodes are often out of the communication range of 
the receiver, therefore the channel detection failure may not lead to collisions. 

C. Multi-hop Topologies with Multiple End-devices 

We consider the multi-hop topology in Fig. Q]:). We use the same MAC and physical layer 
parameters as in the single-hop case. We consider the end-to-end reliability as the routing metric 
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and study the performance of the network as a function of the traffic Aj = A, i — 1, N, in the 
range 0.1 — 10 pkt/s, the spread of the shadow fading in the range a = — 6. Moreover, we show 
results for different values of the Nakagami parameter k = 1 — 3 and threshold b = 6, 10, 14 dB. 

In Fig. \T7\ we report the average end-to-end reliability over all the end-devices by varying the 
node traffic rate. The results are shown for different values of Nakagami parameter k with the 
shadowing spread set to a = 6. The impact of the Nakagami parameter k seems more prominent 
than variation of the traffic. Fig. [T8] shows the end-to-end reliability by varying the spread a for 
different values of b. Differently with respect to the single-hop and linear topologies, a variation 
of the outage threshold b has a strong impact on the reliability also for small to moderate 
shadowing spread. In fact, due to the variable distance between each source-destination pair, the 
fading and the outage probabilities affect the network noticeably. Nonetheless, this effect is well 
predicted by the analytical model we have developed. 

VII. Conclusions 

In this paper, we proposed an integrated cross-layer model of the MAC and physical layers for 
unslotted IEEE 802.15.4 networks, by considering explicit effects of multi-path shadow fading 
channels and presence of interferers. We studied the impact of fading statistics on the MAC 
performance in terms of reliability, delay, and power consumption, by varying traffic rates, inter- 
nodes distances, carrier sensing range, and SINR threshold. We observed that the severity of the 
fading and the physical layer thresholds have significant and complex effects on all performance 
indicator at MAC layer, and the effects are well predicted by the new model. In particular, the 
fading has a relevant negative impact on the reliability. The effect is more evident as traffic and 
distance between nodes increase. However, depending on the carrier sensing and SINR thresholds, 
a fading with small spread can improve the reliability with respect to the ideal case. The delay 
for successfully received packets and the power consumption are instead positively affected by 
the fading and the performance can be optimized by opportunely tuning the thresholds. 

We believe that the design of future WSN-based systems can greatly benefit from the results 
presented in this paper. As a future work, a tradeoff between reliability, delay, and power 
consumption can be exploited by proper tuning of routing, MAC, and physical layer parameters. 
Various routing metrics can be analyzed, and the model extended to the presence of multiple 
sinks. 
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Fig. 1. Example of topologies: single-hop star topology (on the left), multi-hop linear topology (in the center) and multi-hop 
topology with multiple end-devices (on the right). 
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Fig. 2. Reliability vs. traffic rate A for the star network in Fig. QJ) with iV 



= 7 nodes, r = 1 m, a = -76 dBm, 6 = 6 dB. 
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Fig. 3. Delay vs. traffic rate A for the star network in Fig.QJ) with N = 7 nodes, r = 1 m, a — —76 dBm, b — 6 dB. 




Fig. 4. Power consumption vs. traffic rate A for the star network in Fig. [T^) with N 
b = 6 dB. 
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Fig. 6. Delay vs. radius r for the star network in Fig. QJ) with N — 7 nodes, A = 10 pkt/s, a = 



-76 dBm, b = 6 dB. 
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Fig. 9. Delay vs. a for the star network in Fig. [T^) with N = 7 nodes, r = 5 m, A = 10 pkt/s, b = 6 dB. 




Fig. 10. Power consumption vs. a for the star network in Fig. QJ) with N — 7 nodes, r = 1 m, A = 10 pkt/s, b — 6 dB. 
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Fig. 11. Reliability vs. <r for the star network in Fig. QJ) with N = 7 nodes, r = 1 m, A = 10 pkt/s, a — —76 dB. 
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Fig. 12. Reliability vs. a for the star network in Fig. [TJi) with N = 7 nodes, r = 1 m, A = 5 pkt/s, a = —56 dB, 6 = 6 dB. 
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Fig. 13. End-to-end reliability vs. number of hops for the linear topology in Fig. |TJ)) with N = 5 nodes, r = 1 m, A = 2 pkt/s, 
a = -76 dB, b = 6 dB. 




Fig. 14. End-to-end reliability vs. hop distance r for the linear topology in Fig. [Tp) with iV 
a = -76 dB, b = 6 dB. 
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Fig. 15. End-to-end reliability vs. a for the linear topology in Fig.[Tp) with N = 5 nodes, r = 1 m, A = 2 pkt/s, 6 = 6 dB. 
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Fig. 17. End-to-end reliability vs. traffic rate A for the multi-hop topology in Fig. [TJ;) with a — —76 dB, b — 6 dB, a = 6. 




Fig. 18. End-to-end reliability vs. <r for the multi-hop topology in Fig. QJ) with A = 2 pkt/s, a — 



-76 dB. 



